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Xspect, estimation of the angular power spectrum by 
computing cross-power spectra with analytical error bars 
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ABSTRACT 

We present Xspect, a method to obtain estimates of the angular power spectrum of the 
Cosmic Microwave Background (CMB) temperature anisotropics including analytical 
error bars developed for the Archeops experiment. Cross-power spectra are computed 
from a set of maps and each of them is in itself an unbiased estimate of the power 
spectrum as long as the detector noises are uncorrelated. Then, the cross-power spectra 
are combined into a final temperature power spectrum with error bars analytically 
derived from the cross-correlation matrix. 

This method presents three main useful properties : (1) no estimation of the noise 
power spectrum is needed, (2) complex weighting schemes including sky covering and 
map noise properties can be easily taken into account, and corrected for, for each input 
map, (3) error bars are quickly computed analytically from the data themselves with 
no Monte-Carlo simulations involved. Xspect also permits the study of common fluc- 
tuations between maps from different sky surveys such as CMB, Sunyaev-Zel'dovich 
effect or mass fluctuations from weak lensing observations. 
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1 INTRODUCTION 

The measurement of the angular power spectrum 
of the CMB anisotropics, GiS, has become one of 
the most important tools in modern cosmology. 
As long as they remain in the linear regime, the 
fluctuations p re dicted by most inflationary scenari i 
<Hu et al. 19971: iLinde et ah 19991: iLiddle fc Lvth 200d) 
lead to Gaussian anisotropics on the CMB. Thus the 
angular power spectra in temperature and polarization 
contain all the cosmological information on the CMB sky. 
Cosmological parameters and other physical quantities of 
interest in the early Universe can be directly derived from 
them. In parallel to the ex plosion of CMB dat asets both in 
size and quality (W MAP teennett et al. 2003[). Archeop s 
jBenoit et al. 20031). Bo omeran g fctuhl et al. 20o3) . 
Maxim a (Ijaffe et al. 2003T) DASI faalverson et al. 200S ). 
VSA iGrainge et al. 20031) . CBI dSievers et al. 2003) . 
ACB AR lll^ioet^d^OOMl) ) . fast codes have been developed 
to esti mate the CMB angular power spectrum (C MB- 
FAST iSeliak fc Zaldarriaga 19961 IZaldarriaea et al. 199st 
IZaldarriaea fe Seliak 2000T) CAMB jLewis et al. 200d) ) 
allowing us to compare fast and efficiently theory 
and obs ervations us ing powerful statistical tests (CM - 
BEASY JPoran 20031) . COSMOMC jLewis fc Bridle 20021) . 
JPouspis et al. 200 J) ). Furthermore, huge efforts are under- 



taken to ease the estimation of the angular power spectrum 
from input CMB maps in order to cope with larger, deeper 
and more complex sky surveys in a reasonable amount of 
computing time. 

Excluding very specific methods - for example those 
which are under study for the Planck satellite mission 
and which take advantage of the P lanck ring scanning 
strategy (Ivan Leeuwen et al. 20021 : IChallinor et al. 20021 : 
lAnsari et al. 20031) - most CMB power spectrum estima- 
tors can be grouped into two categories : maximum like- 
lihood and 'pseudo'-Ci estimators. A complete review and 
comparison bet ween the two methods can be found in 
lEfstathiou 2004 here we just discuss the key points of each 
of them. 

Maximum likeli hood methods for tempera- 
ture anisotropies feond et al. 199^ : iTeemark 19971 : 



Borrill 1999a) are based on the maximization of the 
quadratic likelihood. The method estimates the sky angular 
power spectrum from the angular correlation function of 
the data. Error bars for the power spectrum are generally 
computed directly from the likelihood function which is 
either fully sampled in the range of interest or approxi- 
mated by a quadratic form. Dealing with an inhomogeneous 
coverage of the sky involves a great computational com- 
plexity 0(Np ix ^, where N P i X is the number of pixels of 
the input map. Therefore, these methods are very CPU 
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time consuming for current lar ge datasets li ke WMAP 
and probably not well adapted (Bor rill 1999bl) for future 
satellite missions like Planck which will produce maps of 
the sky of more than N P i x = 5 x 10 7 pixels. A generalization 
of these me thods to the analysis of CMB pol arization is 
discussed in lTeemark fc de Oliveira-Costa 200ll 

Alternatively, the so-called 'pseudo-CV 's estimators 
compute directly the 'pseudo ' angular power spectrum from 
the data. Then, they correct it for the sky coverage, beam 
smoothing, data filtering, pixel weighting and noise biases. A 
comprehensiv e description of this method was first given by 
IPeebles 19731 and an applica tion to the angular clu stering 
of galaxies can be found in IPeebles fc Hauser 1974 More 
recently, several approaches to this method have b een de- 
veloped. Among them, S PICE ij Szc^aud i et al. 200 lT) and its 
extension to polarization jChon et al. 2004) compute in the 
real space first the correlation function to correct for the sky 
coverage bias and then the power spectrum from the latter. 
A pseudo-CV estimator in the spherical harmonic space ap- 
plie d to CMB experim ents is given in IWandelt et al. 200 ll 
and lHivon et al. 20021 (MASTER). They computes directly 
the power spectrum before correct it for the different biases. 
An appro ach applied to apod ised regions of the sky is pre- 
se nted in lHanse n et al. 20021 and extended to polarization 
in lHansen fc Go rski 2003J. These estimators can be evalu- 
ated using fast spherical harmonic transforms O {n^(^\ and 

therefore provide fast and accurate estimates of the CVs. 
However, they require an accurate knowledge of the instru- 
mental setup and noise in order to correct them for the biases 
discussed previously. In fact, they use an estimation of the 
power spectrum of the noise in the map, generally computed 
via Monte-Carlo simulations, which is subtracted from the 
original power spectrum. This is also used to estimate the 
error bars in the power spectrum by calculating the variance 
of the Cis over the set of simulations. 

In this paper, we describe a method to estimate the Ces 
by computing the cross-power spectra between a collection 
of input maps coming either from multiple detectors of 
the same experiment or from different instruments. The 
'pseudo' cross-power spectra are explicitly corrected for 
incomplete sky coverage, beam smoothing, filtering and 
pixelization. Assuming no correlation between the noise 
contribution from two different maps, each of the cor- 
rected cross-power spectra is an unbiased estimate of the 
Ces. Analytical error bars are derived for each of them. 
The cross-power spectra, that do not include the classical 
aitfo-power spectra, are then combined using a Gaussian ap- 
proximation of the likelihood function. In the same way, we 
can also compute the estimate of the common angular power 
spectrum, (j^ ommon ( f s ky m aps from different experiments. 

A similar method also based on the combination of a set 
of cross-power spectra has been first u sed to obtain recent re - 
sults from the first year WMAP data Jffinshaw et al. 2003). 
The main difference between the method presented is this 
paper and the WMAP one is the determination of the cross- 
correlation matrix (see Sect[3J of the corrected cross-power 
spectra used both for the combination of these into a single 
power spectra and for the estimation of the error bars on 
the latter. The WMAP team estimates the cross-correlation 
matrix from a model of the data. This includes specific terms 



related to the WMAP data such as the contribution from 
point sources and the uncertainties on the beam window 
functions as well as a term related to the CMB anisotropies 
which is estimated from a fiducial model. The WMAP cross- 
correlation matrix, used for the combination of the cross- 
power spectra, does not incorporate the effects of mode cou- 
pling. In a further step, they account for the mode coupling 
and the dependence on the fiducial CMB model for the com- 
putation of the uncertainties on the final power spectrum. 

By contrast, the method presented here computes the 
cross-correlation matrix directly from the cross-power spec- 
tra estimated from the data. This allows us to include nat- 
urally the mode coupling in this matrix. Further, this per- 
mits the computation of analytical error bars (as described 
in Sect. |3J which are very compatible with those obtained 
from simulations (see Sect. Because of the above, this 
method can be applied without modification to the estima- 
tion of the power spectrum of the correlated signal between 
a set of maps of the sky coming from multiple instruments 
with potentially different sky coverages. For example, we 
have used this method on Archeops data for the estima- 
tion of the CMB angular power spe ctrum and the contr i- 
bution from foregrounds to this one jTristram et al. 20051) . 
We have also used it for the estimation of the foreground 
emission at the sub-millimeter and millimeter wavelength by 
cross-correlating the Archeops data with foreground dust 
templates. 

In Sect.|5J we remind to the reader the computation of 
the cross-power spectra from 'pseudo' cross-power spectra. 
In Sect. |3 we specify the correlation between cross-power 
spectra and between multipoles. Analytical expressions for 
the error bars and the covariance matrix for each cross- 
power spectra are derived. Section 0] discusses the combi- 
nation of the cross-spectra from either a single full data set 
(Sect. rOl or several independent experiments (Sect. FOll . 
Finally, Xspect is applied to simulations of the Archeops 
balloon-borne experiment in Sect.|3 



2 CROSS-POWER SPECTRA 

Under the assumption of uncorrelated noise between detec- 
tors, the cross-power spectrum computed from sky maps 
is a non noise-biased estimate of the angular power spec- 
trum. In general, when computing the Ces, other instrumen- 
tal effects as beam smoothing, incomplete sky coverage and 
time ordered data filtering need to be taken into account. 
These effects and the way to correct them have been deeply 
covered in the literature for classical 'pseudo-CV estimators 
JSzapudi et al. 200ltlHivon et al. 20021: lHansen et al. 20021) 
which are noise biased as they use directly the auto power 
spectrum of each detector map. As shown in the following, 
these corrections can be extended to the case of cross-power 
spectra. 

The CMB temperature anisotropies, AT, over the full- 
sky can be decomposed into spherical harmonics as follows, 

Ar(n)=^a, m ri m (n) (1) 

£m 

where the coefficient at m are given by 

ai m = / AT(h)Y em (n)(Kl . (2) 
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The CMB temperature field AT 



redicted by most 
" Linde et al. 19991: 



inflationary mo dels dLiddle fc Lvth 200(1 [ 
iHu et aL 1 997'l is in general Gaussian distributed so that 
the ensemble average of the ae. m coefficients are 



(a em ) = (3) 
(ae m ae'm'} = (Ct)8 U i8 mm i (4) 

An unbiased estimate of the CMB temperature power 
spectrum (Ci) is therefore given by 



C, 



k'tti 



E 2f -;- !. 



(•») 



Ground-based and balloon-borne CMB experiments 
present an inhomogeneous sky coverage. On the other hand, 
satellite experiments like COBE, WMAP and Planck, al- 
though they provide full-sky maps, residuals of foreground 
contamination in the Galactic plane and point sources con- 
tamination make impossible to use the complete maps when 
computing the CMB angular power spectrum. Furthermore, 
for most CMB experiments, the noise properties vary con- 
siderably due to different redundancies between pixels of the 
same map. So obtaining an estimate of the power spectrum 
requires a differential weighting of pixels within the same 
map which translates into an effective inhomogeneous sky 
coverage. 

The decomposition in spherical harmonics of the ob- 
served temperature anisotropics including weighting can be 
written for a single detector as follows 



di 



dnr™ p (h)w(h)Y em (n) 



N PIX 

J2^pT^ p W(n p )Y em (n p ), 
P =i 



(6) 



(T) 



where Q~ 



for equal pixels (as in HEALPix pixeliza- 



tion, iGorski et al. 1993) and W is the mask applied to the 
input sky temperature T map . This temperature can be de- 
composed into signal and noise, T map = T st9nal + T" oise 
which are assumed to be uncorrelated, i.e. 



(T 



• ignal rpnoise 



(8) 



The effect of a non-homogeneous co verage of the sky 
can be described in spherical harmonics (IPeebles 1973T) by 
a mode-mode coupling matrix Mui which depends only on 
the angular power spectrum of the weighting scheme (here- 
after weighting mask) applied to the sky to account for the 
incomplete coverage, the removal of the Galactic plane and 
the inhomogeneous noise properties of the detector. 

Thus, the 'pseudo-CV estimator Di is defined as follows 



D e = 



2! 



— y 

+i ^ 



(9) 



The relation that links the 'pseudo' power spectrum, 
directly measured on the sky, and the power spectrum Ce of 
the CMB anisotropies is given by 



D t =22 M a'\Pt' B t'\ 2F t'( C e') + W ( 10 ) 
t' 

where Be is the beam transfer function describing the beam 



smoothing effect; pi is the transfer function of the pixeliza- 
tion scheme of the map describing the effect of smoothing 
due to the finite pixel size and geometry; Ft is an effective 
function that represents any filtering applied to the time 
ordered data; and (Ne) is the noise power spectrum. 



2.1 Pseudo cross-power spectrum 

An unbiased estimate of the cross-power spectrum Cf B be- 
tween the full sky maps of two independent and perfect de- 
tectors A and B can be obtained from 



r<AB 



21 + 1 



E 



B « 



(11) 



where af m and af m are the coefficients of the spherical har- 
monic decomposition of maps A and B respectively. 

In the same way, we can compute the 'pseudo' cross- 
power spectrum Df B between any two detectors A and B 
by generalizing Eq. 1101 



Di 



\ 2 Bf;Bf,F£ B (cr) + {Nt B ) . (12) 



Each of the terms in Eq. 1121 is described in more details in 
the following subsections. 



2.2 



Noise cross-power spectrum, Nf B 



The main advantage of using cross-power spectra is that the 
noise is generally uncorrelated between different detectors 



/ A B 
\ n lm n t'm' 



(13) 



This assumption will be maintained throughout this paper. 

Thus, the cross-power spectra are straightforward esti- 
mates of the angular power spectrum on the sky and for two 
different detectors (i.e. A ^ B), the 'pseudo' cross-power 
spectrum reads 



D AB = 



\cP) 



(14) 



2.3 Coupling kernel matrix, M U B 

The coupling kernel matrix Mm/, introd uced in Eq. 1101 and 
described in details in lHivon et al. 20021 reads 



M «' = ^ i E( 2£ " + 1 ) w ^ 



£ £' 








(15) 



where We = 2 i+i ^2 m =-i \ We ™\ 2 ^ s l ne P ower spectrum of 
the mask. It takes into account the mask applied to the 
data where mask represents both the sky coverage and the 
weighting scheme. 

Equation ll5l can be easily generalized to the case of two 
different masks applied respectively to each map of the two 
detectors involved in the cross-power spectrum calculation. 
Replacing the quadratic terms \we m \ 2 by (w^wf^j in the 
computation of Eq. I15l leads to 



Ml, 



2t! + 1 

47T 



^(2£" + l)Wt, 



I" 




(16) 
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where Wf B = jzpL Y] m __ e wf m wf^, the cross-power spec- 
trum of the masks. 

This property allows us to deal with independent masks 
representing different sky coverages and to apply an appro- 
priate specific weighting scheme to each detector map. Note 
that the correction in the multipole space discussed here is 
fully equivalent to an appropriate normalization of the cross 
correlation between the two sky masked maps in real space. 
This analogy is important as it helps to understand why no 
fully overlapping masks for the input maps can be consid- 
ered. 



2.4 Filter function, Ff 

The filter function Fi accounts for the filtering of the time 
ordered data which is generally needed in most CMB exper- 
iments either to avoid systematic effects or to reduce cor- 
related low frequency noise. The time domain filtering is 
performed along a preferred direction on the sky (scanning 
direction) and so leads commonly to an anisotropic sky even 
if the assumption of initial isotropic temperature fluctua- 
tions holds. In this case, the estimates of the angular power 
spectrum provided by Eg. HUI and Eg. 1141 are not exact any 
more and should be corrected for a function both in I and 
m, Fe t7n . Obtaining accurate estimates of such a correction 
is particularly com plex and for most cases, as proposed by 
iHivon et al. 2002 . the correction for an effective Fe is good 
enough for the accuracy required in the reconstruction of the 
CMB power spectrum. The Ft function can be, for exam- 
ple, computed via Monte-Carlo simulations of the sky from 
which mock time ordered data are produced for each of the 
detectors involved and then filtered. 

From an initial theoretical CMB power spectrum, we 
compute a large number of realizations of the sky using the 
HEALPix software synfast iGorski et al. 199Sl) and compute 
mock time ordered data from the scanning strategy of each 
detector. Maps are then computed with and without filter- 
ing before re-projection. The Ft function is obtained from 
the mean ratio of the 'pseudo' power spectra of the filtered 
and not filtered maps. The latter are obtained using the 
HEALPix software anafast. 

In the case of the cross-power spectra and consider- 
ing the previous approximation, an effective filter function 
Ff B — \fpfFf will be considered in the following. As de- 
fined, the effective filtering function allows us to consider 
detectors for which the time domain filtering is different. 

Note that, for nearly white noise and all-sky surveys 
such as WMAP or Planck missions, filtering may not be 
required and thus Ft — 1. 



JSouradeep fe Ratra 200ll : lFosalba et al. 20021) . The work of 
IWandelt fe Gorski 200 J presents how to convolve exactly 
two band limited but otherwise arbitrary functions on the 
sphere - which can be the 4-7T beam pattern and the sky 
descriptions. An analytic framework for studying the lead- 
ing order effects of a non-circular b eam on the CMB power 
spectrum estimation is proposed in lMitra et al. 20041 

The authors of this paper present Asymfast 
jTristram et al. 20041) . a general method to estimate 
an effective Bg function taking into account the asymmetry 
of the main beam and the scanning strategy. Asymfast 
is based on a decomposition of the main beam pattern 
into a linear combination of Gaussians which permits fast 
convolution in the spherical harmonic space along the 
scanning strategy. 

2.6 Cross-power spectrum from 'pseudo' 
cross-power spectrum 

The cross-power spectrum Cf B can be obtained from the 
'pseudo' cross-power spectrum D AB by resolving Eq. 1121 
which leads to invert the coupling kernel matrix Mu' AB . 
In general for complex sky coverage and weighting schemes, 
this matrix is singu lar and can not be inverted directly. To 
avoid this problem, IHivon et al. 20021 proposed the binning 
of the coupling kernel matrix which reduces considerably 
the complex correlation pattern in the Ces introduced by 
the applied mask. The binning is obtained by applying the 
operators Pu and Qn, as follows 



c b AB 


= PuCf B 


(17) 


dJ b 


= P w Df B 


(18) 


Qlb 


- p- 1 


(19) 



The solution of Eo. 1141 in the new base reads 
Cp'=M£Dp> (20) 



with 

Mw = Pu ( M AB F(i Pfi B A B B ) Q llv . (21) 

As the correlation between multipoles I depends very 
much on the instrumental setup, the sky coverage and the 
weighting scheme, the binning has to be defined for each 
experiment. It is a compromise between a good multipole 
sampling and low correlations between adjacent bins. 

Hereafter to avoid confusion and makes the notation 
simpler all equations will be written in £ instead of b. 



2.5 Beam window function, Bf 

The beam window function Bt describes the smoothing ef- 
fect of the main instrumental beam under the hypothesis 
of circularity. The latter does not hold in general as for 
most experiments the main beam pattern is asymmetric. 
As beam uncertainties have become the most important 
source of systematic errors, taking into account the asymme- 
try of the beam pattern is necessary. Sever al solutions have 
been proposed ei ther circularizing the beam l|Wu et al. 2000t 
iPaee et al. 20031) or assuming an elliptical Gaussian beam 



3 CROSS-CORRELATION MATRIX - 
ANALYTICAL ERROR BARS AND 
COVARIANCE MATRIX 

From N input maps we can obtain N(N — l)/2 cross-power 
spectra C AB (A 7^ B) which are unbiased estimates of the 
angular power spectrum but which are obviously not inde- 
pendent. In this section, we describe the estimation of the 
cross-correlation matrix between cross-spectra and between 
multipoles. We show how the error bars and the covariance 
matrix in multipole space can be deduced for each cross- 
power spectra. 



Xspect 5 



3.1 The cross-correlation matrix, S 

Given a sky map from a detector A, it can be combined to 
each of the other detector maps to form N — 1 cross-power 
spectra which will therefore be highly correlated. Further- 
more, due to the masking we also expect that each cross- 
power spectra will be correlated for adjacent multipoles and 
thus, correlations between adjacent multipoles will be also 
present between different cross-power spectra. To describe 
this complexity we define the cross-correlation matrix 



~AB,CD 

1 U> 



={(ct B -{ct B )){c$ D -{c?, D )y) 



of the cross-power spectra AB (A 7^ B) and CD (C 7^ D) 
which can be fully computed as shown in the following. 

From Eq. D~D we can express the 'pseudo' cross-power 
spectrum between detectors A and B as follows 



D. 



M? e B Cfi B 



(22) 



where M AB = Mf^EfiEf, and E e = 



Bty/Fi and there- 
fore the corrected cross-power spectrum for detectors A and 
B is given by 



Cf B = (M% , y 1 D$ B . 



(23) 



Using the above expression and following the previous 
definition, the cross-correlation matrix H..f ' CD reads 



■i B > CD EE (AC? B AC?, D ' 



= (Mfi- 1 *Dt?(M%°- i AI% D )*) 

= M^-^AD^ADg^^M^-Y (24) 

The above equation as it stands can not be used in 
practice for £ above ~ 10 as the calculation of 



r AC r BD K A 7^0 K 



mm' Z\m\ ^2 m 2 



(25) 

described in Appendix lAl is numerically unstable 
flVarshalovich et al. 198sIl However, for high multipoles and 
sufficiently large sky coverage (as it is the case for satellite 
missions as WM AP and Planck) it can be simplified 
l)Efstathiou 2004) by replacing C AC and C BD by Cp c 
and C BD respectively and then applying the completenes s 
relation for spherical harmonics dVarshalovich et al. 198 j) . 
This is because the Ke m ti m i, which is diagonal for full sky 
coverage, is quasi-diagonal for large sky coverage. 

Under the above hypothesis and following Appendix lAl 
the cross-correlation matrix reads 



M 



AC,BD\ r<AC nBD 



2£ 2 + l 
2£ 2 + l 



where 

M% 2 (W AC - BD ) = E&EZm™ (w ac ' bd ) Eg El 
and is the quadratic coupling kernel matrix 



(212 + 1) 

47T 



AB,CD\ 



^2(213 + 1)W£ 



13 



,cd( 11 12 13 




(27) 



associated to the cross-power spectrum of the product 
of the masks W AB ' CD = ^E,,'""'^™'. The 
W em represents the spherical harmonic coefficients for 
the product of the masks associated to detectors A and B. 

Equation 1261 can be further simplified by assuming uni- 
form weighting and the same sky coverage for all detectors, 
as well as a diagonal dominated coupling kernel matrix, 



^AB,CD 



±\c AC c BD + c AD c BC ] 



(28) 



In this case, the effect of a non-homogeneous sky cov- 
erage is represented by a simple function vi which can 
be associated to the effective number of degrees of free- 
dom in the \\„ d istribution of the Ces over the sky (see 
iHivon et al. 2002ft 



(2£ + 1)A, 



W4 



where Wi is the i-th moment of the mask 



4tt 



<mw\n) . 



(29) 



(30) 



In Eo. l26l and l28l the Ci can be either cross-power spec- 
tra or auto-power spectra depending on the combination of 
A, B, C and D (the only condition is A 7^ B and C 7^ D). In 
one hand, noise terms, as included in the auto-power spec- 
tra, appear in the analytical form of the correlation matrix 
S. On the other hand, for a set of 4 independant detectors 
(A 7^ B 7^ C 7^ D), the correlation matrix is the variance of 
the signal. 



3.2 Covariance Matrix and error bars associated 
to the cross-power spectra 

It is important to notice that the cross-correlation matrix 
contains all the needed information to perform error bars 
and covariance matrix for each single cross-power spectrum. 
From Eg. 1261 (or Eg. 1281 for the approximated form), the co- 
variance matrix providing the correlation between adjacent 
multipoles is given by 

-AB,AB 



Gov AB (£,£') = 



1 



\cf B C AB +C AA C BB ] 



(31) 



In the same way, we can also write the variance of 
each cross-power spectra which corresponds to the error bars 
ACe AB associated to each cross-power spectra 



(AC e AB ) 2 



1 



AB\ 2 . ^AA^BB 
) + Of Of 



(32) 



In these cases, instrumental noise appears in the auto- 
power spectra Cf x (for which N xx 7^ 0). As it is computed 
directly from the data, no estimation of the noise is needed. 
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4 COMBINED ANGULAR POWER 
SPECTRUM 

This section describes the estimation of the final CMB an- 
gular power spectrum and the error bars associated to it 
using the corrected cross-power spectra described in Sect. [5] 
in addition to the cross-correlation matrix computed in 
Sect. [3] We propose a simple but efficient way to com- 
bine them by maximizing a quadratic likelihood function 
and to deduce the final error bars from the cross-correlation 
matrix. For this paper we have considered a very general 
form for the likelihood function limited to a simple CMB 
plus noise model. However, the likelihood function can be 
adapted more specifically to the data, for example to in- 
clude point sources as it is the c ase for the WMAP results 
presented in Hinshaw ct al. 2003. More generally, it could be 
exte nded to algorithms of sp ectral matching decomposition 
{e.g. iDelabrouille et al. 2 003). In the following, we present 
how this method can be used to provide both an estimate 
of the Ces from a full data set (Sect. lCT and of the angular 
power spectrum of common sky fluctuations between maps 
coming from two or more independent surveys (Sect. [Ot . 

4.1 Gaussian approximated linear combination of 
the cross-power spectra 

Once all possible cross-power spectra have been computed 
from a data set, we dispose of N(N — l)/2 different but 
not independent measurements of the angular power spec- 
trum, C'is. To combine them and obtain the best estimate 
of the power spectrum Ct, we maximize the Gaussian ap- 
proximated likelihood function 



auto-power spectra (including noise). The variance of final 



21n£ : 



[(CI 



Ce)\E \y e ,(C 3 e , - C t i) 



(33) 



where |E«/| 1 -' = E,f e f' CD is the cross-correlation matrix 
of the cross-power spectra described before ( i and j £ 
{AB,A 7^ £>}). The auto-power spectra are not considered. 

From this and neglecting the correlation between adja- 
cent multipoles, it is straightforward to show that the esti- 
mate of the angular power spectrum is 



lEjIg-^ + QIs- 1 ^] 
2 EJs-MJi 



(34) 



The final covariance matrix can be obtained from 
Eg. 1331 



Cov(£,£') 



V lE- 1 !^, 



and the final error bars are given by 
ACe) = ' 



(AC/] 



-i i»j 



(35) 



(36) 



Depending on the cross-correlation between cross-power 
spectra, the instrumental noise variance is reduced by a fac- 
tor comprised between N, the number of independant detec- 
tors, and N(N — f)/2, the number of cross-power spectra. 

For the noise dominated case, the correlation between 
different cross-power spectra can be neglected and the cross- 
correlation matrix S becomes diagonal (see eg . 1261 and I28B . 
The values in the diagonal are the 7V(iV — l)/2 products of 2 



power spectrum is then proportional to 



In any 



JV(iV-l)/2 ' 

case, when combining cross-power spectra, the upper limit 
for the variance comes from the combination of N indepen- 
dent detectors (proportional to -jj). 



4.2 Common angular power spectrum 

The Xspect formalism allows us to compare two or more 
different sets of sky maps coming from two or more 
independent experiments. In this respect, the quantity of 
interest is the power spectrum of the common fluctuations 
on the sky with the same physical origin. We will call 
the latter common angular power spectrum (J1° mmon for 
simplicity. In fact, if we compare for example two sets of 
CMB maps at low and high frequencies, the foreground 
contamination will be different and we can expect to obtain 
a better estimation of the CMB power spectrum. In the 
same way, for two different experiments, systematic effects 
will be decorrelated and will not contribute to the common 
angular power spectrum calculated from them. In addition, 
template maps of foreground emission can be correlated to 
the CMB experiments maps to monitor and to subtract 
foreground residuals. 

To simplify the notation we will only consider two sets 
of sky maps A and B, corresponding to Na and Nb detec- 
tors respectively. Following the previous paragraph consider- 
ations, by cross-correlating each detector map from A which 
each detector map from B and correcting the 'pseudo' cross- 
power spectra as described before, we can form Na X Nb 
cross-power spectra which are unbiased estimates of the 
common power spectrum but not fully independent. 

The cross-power spectra formed in this way can be 
noted C e % ' 3 and their cross-correlation matrix H^, 3 ' " " 
reads 



M 



A-iBj -1 



A U B V 



A u Bv -1\T 



(37) 



which can be computed using previous approximations. 

Finally, eo uations I34l and 1361 are used to obtain the esti- 

-common 

mate of the common angular power spectrum Ct and 
the error bars associated to it. 

This method allows us to deal naturally with very dif- 
ferent experimental configurations such as for example com- 
pletely different scanning strategies, different timeline filter- 
ing and different beam patterns for each of the detectors 
involved. Furthermore, it can be easily generalized to work 
with non CMB data such as maps of the Sunyaev-Zel'dovich 
effect in clusters of Galaxies or maps of mass fluctuations 
from weak lensing observations. 



5 XSPECT APPLIED TO ARCHEOPS 
BALLOON-BORNE EXPERIMENT 

Xspect was mainly developed to measure the CMB tem- 
perature angular power spectrum from the Archeops data 
as well as to compare these data to other observations of 
the sky as for example those from the WMAP satellite. 
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Archeops feenoit et al. 20031) is a balloon-borne experi- 
ment conceived as a precursor of the Planck High Fre- 
quency Instrument 1 . It consists of a 1.5 m off-axis Gregorian 
telescope and of an array of 21 photometers cooled down 
to ~ 100 mK and which operates in 4 frequency bands: 
143 and 217 GHz (CMB channels) and 353 and 545 GHz 
(Galactic channels). For the latest flight campaign, the en- 
tire Archeops data cover about ~ 30% of the sky, including 
the Galactic plane. 



5.1 Simulations of the Archeops data set 

As a first step in the application of Xspect to the Archeops 
data set, we have produced 500 simulations of the sky ob- 
served by Archeops including both the CMB signal and the 
detector noise for the six most sensitive Archeops detectors 
at 143 and 217 GHz. 

These simulations are computed from realizations of the 
CMB sky at nside = 512 (corresponding to a pixel size of 
~ 7 arc min) for the Arch eops best-fit CMB model pre- 
sented in lBenoit et al. 20031 . The sky maps are convolved by 
the main beam pattern of each of the Ndetec = 6 Archeops 
detectors using the beam transfer function and then depro- 
jected following the Archeops scanning strategy to pro- 
duce mock Archeops timelines. The beam transfer func- 
tions were computed individually for each detector from 
Jupiter's cr ossings in the data us ing the Asymfast method 
presented in lTristram et al. 20041 

The detector noise is computed from the time power 
spectrum of each of the detectors and added to the mock 
signal timelines. The method used for the estimation of 
the time noise power spe ctra is described in details in 
lAmblard fc Hamilton 2004 To suppress 1//- noise, remain- 
ing atmospheric and galactic contamination as well as non- 
stationary high frequency noise, the Archeops time streams 
are Fourier filtered using a bandpass filter between 0.1 and 
38 Hz. We apply the same filtering to the simulated time- 
lines before projection. The corresponding multipole filter- 
ing function Fi was independently calculated by specific sim- 
ulations. The filtered timelines are then reprojected using 
the Archeops scanning strategy to produce Ndetec simu- 
lated co-added maps for each simulation run. 

A Galactic ma sk deduced from a SFD IRAS map 
flSchleeel et al. 1 998?). extrapolated to 353 GHz, using a cut 
in amplitude (greater than 0.5 MJy.str) is applied to the 
simulated maps. This reduces the total Archeops coverage 
to ~ 20% of the sky. We have used two different weighting 
schemes on the maps. The first one is an uniform weight- 
ing and the second one, a l/rjp ix weighting per pixel and 
per detector, where Op ix is the noise variance for the given 
pixel. 

Then we apply the Xspect method to each simulations- 
run. The Ndetec(Ndetec — l)/2 cross-power spectra are com- 
puted with their associated cross-correlation matrices and 
combined as described in Sect. |4~T1 using the Gaussian ap- 
proximation of the likelihood function. 

In addition to the above simulations, we produced, in 
the same way, 500 non noisy simulations from which we can 
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Figure 1. Mean angular power spectrum obtained from 500 sim- 
ulations of Archeops observations of the C MB sky (see text for 
details). We use the Archeops best-fit model iBenoit et al. 2003) 
to produce CMB maps which are smeared out using the beam 
pattern of each of the six most sensitive Archeops detectors and 
deprojected into mock Archeops timelines. Noise is then added 
up for each detector and the resulting timelines projected into 
sky maps are analyzed using Xspect. The error bars shown here 
are computed analytically as described in previous sections. The 
Xspect estimate of the angular power spectrum (in red) is an 
unbiased estimate of the input CMB model (in yellow). 

extract for example the sample variance for the Archeops 
coverage. 

Figure shows the mean angular power spectrum com- 
puted from the angular power spectrum estimates obtained 
for each of the 500-simulations runs. The error bars are an- 
alytically computed. We observe that the Gaussian approx- 
imated combination of cross-power spectra is an unbiased 
estimate of the input CMB model for the angular power 
spectrum. This mean angular power spectrum is a combi- 
nation of two weighting schemes. Up to £ — 90, we use 
uniform weighting whereas, for larger £, a mask inversely 
proportional to the noise in each of the pixels has been ap- 
plied. This allows us to choose the smaller error bars in each 
region. 

Figure [5] shows analytic estimates of the error bars in 
the final angular power spectrum compared to the dispersion 
for each multipole bin over the 500 Monte-Carlo simulations. 
The Archeops sky coverage is very inhomogeneous due to 
the particular choice of the scanning strategy which tries to 
maximize the area of the sky observed (~ 30%) in a very 
reduced amount of observation time (less than 12 hours). 
Archeops performs large circles on the sky that leads to a 
quasi ring-like coverage with a large uncovered region at the 
center of the ring as shown on Fig. [3] This implies that the 
approximation used to obtain Eg. 1281 is not valid, especially 
at low multipoles, and therefore the error bars obtained an- 
alytically (red solid line) are underestimated up to ~ 20% 
at very low L with a mean of ~ 10%. By computing Ef ( f' CD 
as given by Eg. 1261 the analytic error bars (yellow solid line) 
increase as expected and fit much better the dispersion in 
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Figure 2. Error bars for the ArcheopS data computed an- 
alytically, from Eq. 1281 (red solid line) and from Eq. 1261 (yel- 
low solid line). They are compared to the standard deviation of 
Nsimu = 500 simulations (black dashed line). Differences essen- 
tially originate from the extremely inhomogeneous ArcheopS sky 
coverage (see text for details). 




Figure 3. Sky coverage and weighting scheme applied to the sky 
map (nside = 512) obtained from the ArcheopS most sensitive 
detector at 143 GHz. We observe a ring-like structure with a large 
uncovered area in the center and highly redundant areas on the 
edges of the ring. 



the simulations (black dashed line). The agreement is then 
within 10% over the full range of multipoles with a mean of 
1.5%. Equivalent results are obtained using non noisy sim- 
ulations. 



5.2 Application to the ArcheopS data 

Xspect has been applied to the data from the last ArcheopS 
flight campaign using the six most sensitive bolometers as in 
the simulations presented above. The results of this analysi s 
is presented in a collaboration paper jTristram et al. 2005). 
Cross-correlation with WMAP maps is also under study in 
order to assess the electromagnetic spectrum of the CMB 
anisotropics from 40 to 217 GHz. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we have presented a method, Xspect, for the 
obtention of the CMB angular power spectrum with analyt- 
ical error bars based on four main steps: 

(i) Given N independent input sky maps from different 
detectors either from a single experiment or from multiple 
ones, we estimate N(N — l)/2 'pseudo' cross-power spectra 
as well as N 'pseudo' auto-p ower spectra using the HEALPix 
package iGorski et al. 199a) . 

(ii) Cross-power spectra are obtained by correcting the 
'pseudo' power spectra for weighting scheme, beam smooth- 
ing and filtering as discussed in Section |2] In the same way, 
we compute auto-power spectra with noise. Xspect can deal 
with different and complex weighting schemes for each of 
the detectors involved. The coupling matrix, beam and fil- 
tering transfer functions are precomputed for each pair of 
detectors. 

(iii) We compute the full cross-correlation matrix, 
S,f e f' CD , between cross-power spectra AB (A ^ B) and CD 
(C 7^ D) (Section l3.1H from which we can extract the covari- 
ance matrix and the error bars for each cross-power spectra 
(Sect [£2Jl . 

(iv) Finally, the corrected cross-power spectra are com- 
bined into a single angular power spectrum using their cross- 
correlation matrices which are assumed to be diagonal in 
multipole space. The Gaussian approximation of the likeli- 
hood function is fully justified in the large multipole range. 
In addition, analytical estimates for the error bars and for 
the covariance matrix are computed (Section^ . The covari- 
ance matrix can be used to check the degree of correlation 
between multipoles. 

As this method estimates the angular power spectrum 
using 'pseudo' cross-power spectra, we can obtain a non 
noise-biased power spectrum avoiding the estimation of the 
noise power spectra which requires heavy Monte-Carlo sim- 
ulations. Equally, this permits the estimation of analytical 
error bars which are compatible to those computed from 
simulations. Furthermore, Xspect allows us to obtain the an- 
gular power spectrum of common sky fluctuations between 
two or more experiments. Associated to other surveys used 
as templates, it can provide estimations of systematic resid- 
uals or astrophysical contaminations. Nevertheless, Xspect 
computes only the common structures between maps, as- 
suming a single physical component. This means that, for 
CMB purposes, foregrounds and systematics must be sub- 
tracted beforehand. 

Xspect has been successfully applied to simulations of 
the ArcheopS experiment which presents an inhomoge- 
neous sky coverage and detectors with unequal sensitivities. 
Even for a balloon-borne experiment that can be compared 
to satellite missions neither in noise level nor in sky cover- 
age, the analytical estimates of the errors computed with 
Xspect are accurate within few percents. This property is 
of great interest when producing Monte Carlo simulations 
in the process of testing and improving the estimation of 
the Cis or when checking the robustness of the data analy- 
sis with respect to the various choices of mask, filtering or 
binning. The application of Xspect to the ArcheopS data 
set will be p ublished by the ArcheopS collaboration shortly 
jTristram et al. 20051) . 
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For the WMAP satellite, which covers the full sky in 
a very homogeneous way, it was also interesting to check 
the analytical estimates of the error bars. Xspect analyti- 
cal errors are equivalent to those provided by the WMAP 
team within 15% with a mean of 1.8%. This agreement is 
satisfactory as our analysis was more basic than the one 
used to derive the first year WMAP results presented in 
Hinshaw et al. 2003L In particular all maps at all multipoles 
were used without any point-source specific treatment and 
only two weighting schemes were applied. 

The extension of Xspect to CMB polarization maps is 
under development. As for the temperature power spectrum, 
the Cf E and Cf B power spectra obtained directly from a 
single set of I, Q and U maps via a 'pseudo' power spectrum 
estimator are noise biased. A 'pseudo' cross-power spectrum 
estimator adapted to polarization can solve this problem by 
using independent sets of I, Q and U maps. 
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APPENDIX A: CORRELATION MATRIX 

In this appendix we describe in details the calculation of the cross-correlation matrix Ef e f' CD under the hypothesis of large 
sky coverage which leads to 



E? e ?' CD = (ACf B AC?, D *) (Al) 
= Mi£- 1 (ADf 1 B AD¥ 2 D *yM?,? 2 - 1 f where AX = X - {X) (A2) 
Let us remind that for Gaussian distributed variables Xi with variance (xiXj) = ah the quadri-variance is given by 

(XiXjX k Xl) = (XiXj) (x k Xi) + (XiX k ) (XjXi) + (x t Xi) {XjX k ) 

= o'ijVki + <Hk"]i + oliOjh (A3) 
Using Eg IA3I the main term of the cross-correlation matrix in Eg. lA2l reads 

mi ni2 

/rl A rl B * \/rl C * rl D \ 4- / rl A rl C * \/rl B * rl D \ 4- / rl A rl D \/rl B * rl C * \ 

E\ u l 1 m 1 a l 1 m 1 /\ ll (2m 2 u (2iii2 / ' \ lI *im 1 "£ 2 m 2 / \ a l 1 mi u « 2 m 2 / ' \ "<i mi u £' 2 m 2 /\ u *im 1 %m 2 / 
(2^1 + 1) (2^2 + 1) 

mi m 2 

= -Oil f(2 + _|_ l)(2/2 + 1) ^ < [(^ilmi ^i2m 2 )( c '!lmi c 'i2m 2 ) + (^ilmi ^2m 2 ) (rf;imi^i2m 2 )] 



mi m 2 

where c^ m are the coefficients of the spherical harmonic decomposition of the masked sky map such that 
dim = ^ ' aiim 1 Ke m ii mi (A4) 

l lmj 

with 

Ki imi i 2m2 = ^""^ u>jflil^ imi (6 l i)y^ 2m2 (flj) (A5) 

i 

- E(-«-w ( ^™+'' )' ,J ( ? 5 ( ( * X ™ ) 

fii is the area of pixel i. Wi represents the mask applied to the sky map and we m the coefficients of the spherical harmonic 
decomposition of the mask 

^Vfiiy £m (0i). (A7) 



Wlr, 



Thus, the cross-correlation matrix can be simply written as follows 

{^■Df^ 3 ADl2 CD ^ = (2/1 + 1) (2Z2 + 1) ^ / [{dnm 1 d^2rn 2 ){dnm 1 dl2m2) + (^ilmi ^Hm 2 }(^Hm 1 ^m 2 )] (A8) 

mi m 2 



Let us develop each of the terms in Eg. IA8I using Ee — ptBt\f~F\ 

(^m^'m') ~ ^ ' ^ ' ( a ^imi% 2 m 2 )-E'!l-£'!2^m£imi^'m'< 2 m 2 
^imi £ 2 m 2 



— ^ ' EiiEnKi m ( irni Ki imie i rn i 



Assuming a large sky coverage fefstathiou 2004^ we can replace Cf^ by Cf Y and then, applying the completeness relation 
for spherical harmonics J\ Y£ m (8i)Yi m (8 j) = ^-S(8i — 6j), we obtain 

{df m dj* m ,) = Cf Y EfEj^2 { W ? W Y) ttiYe m (8i)Y e * m ,(8i) 

i 

= CFEfEjK^Z (A9) 
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Thus, under the above hypothesis Eg . IA8I reads 



(ADffADg D - 



P A P C P S r,fl 
^11^11^12^12 

(211 + l)(2/2 + 1) 



E\ J<r (2)AC K {2)BD r AC r BD K {2) AT> „(2) BC r AD r Bc\ /A-irA 
[ A Ilm 1 I2m 2 A ilm 1 I2m 2 L 'll °i2 + A Umi i2m 2 A ilmi i2m 2 U ll °*2 J V Ai Uj 



The above expression is equivalent to that for the coupling matrix M«> presented in lHivon et al. 20021 and can be simplified 
as follows 



AC BD (2^+l)(2f + 1) 



Ers(2)AC T A2)BD \ " \ " \ ~> AC 

mm' mm' {L\mi &2 rn 2 

£ £' h \( £ I' e 2 

oo oy^oo o 

= (2t+l)M$ (W AC - BD ) 
where M' 2 - 1 is the quadratic coupling kernel matrix 



4-7T 



((2£i + 1)(24 + 1)) ' (-1) 



X l 2 ( — 1 \ m l+ m 2 



I I' ll 

m —m' —mi 



rn — jtl 



4 

-rt%2 



(2£ 2 + 1) £ 



(2£ 3 + 1) AB.CD ( il (-2 h 

4tt ^ I 



(All) 



(A12) 



associated to the cross-power spectrum of the products of the masks w* Y = wf * wj for each of the cross-power spectra AB 
et CD, 



W 1 



2f 



—y 

+1 ^ 



AB CD* 



Replacing IaTTI in IaToI we obtain the following expression 



(AD AB AD? 2 D *) 



tjiA t-tC rpij rp-L 



<i (W AC - BD ) Ct?C?°+M<?X (W AD ' BC ) C AD Cg c 



2£ 2 + l 

and thus, using the same abreviation M ( t f l2 (W AC ' BD ) = E^E^M^ (W 



(A13) 



a (?c S/f (2) ( u/ ac,bd\ E g E D tlle cross _ corre l a tion matrix reads, 



M 



Ui 



2£ 2 + l 



(A14) 



